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ABSTRACT 

We report the discovery of 12 new fossil groups of galaxies, systems dominated by a single giant 
elliptical galaxy and cluster-scale gravitational potential, but lacking the population of bright galaxies 
typically seen in galaxy clusters. These fossil groups (FGs), selected from the maxBCG optical cluster 
catalog, were detected in snapshot observations with the Chandra X-ray Observatory. We detail the 
highly successful selection method, with an 80% success rate in identifying 12 FGs from our target 
sample of 15 candidates. For 11 of the systems, we determine the X-ray luminosity, temperature, and 
hydrostatic mass, which do not deviate significantly from expectations for normal systems, spanning 
a range typical of rich groups and poor clusters of galaxies. A small number of detected FGs are 
morphologically irregular, possibly due to past mergers, interaction of the intra-group medium (IGM) 
with a central AGN, or superposition of multiple massive halos. Two-thirds of the X-ray-detected 
FGs exhibit X-ray emission associated with the central BCG, although we are unable to distinguish 
between AGN and extended thermal galaxy emission using the current data. This sample, a large 
increase in the number of known FGs, will be invaluable for future planned observations to determine 
FG temperature, gas density, metal abundance, and mass distributions, and to compare to normal 
(non-fossil) systems. Finally, the presence of a population of galaxy-poor systems may bias mass 
function determinations that measure richness from galaxy counts. When used to constrain power 
spectrum normalization and Sl m , these biased mass functions may in turn bias these results. 
Subject headings: galaxies: clusters: general, galaxies: groups: general, galaxies: clusters: intracluster 
medium, X-rays: galaxies: clusters, surveys 



1. INTRODUCTION 

Fossil groups (FGs) are systems dominated by a sin- 
gle, giant elliptical galaxy, yet their X-ray emission indi- 
cates a deeper cluster-scale gravitational potential. They 
are generally defined as systems with a AR = 2 magni- 
tude difference between the first and second rank galax- 
ies within 0.5r2oo 12 > and they have an extended ther- 
mal X-ray halo with Lx,boi > 10 42 h^Q erg s _1 (Jones 
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et al. 2003). FGs are thought to be old, isolated galaxy 
groups and clusters in which the large galaxies have co- 
alesced through dynamical friction (Ponman et al. 1994; 
Mulchaey & Zabludoff 1999; Jones et al. 2003). This co- 
alesced cluster scenario is further supported by high X- 
ray temperature measurements (up to ~ 4 keV) and by 
the galaxy velocity dispersions (e.g., Mendes de Oliveira 
et al. 2006, 2009; Cypriano et al. 2006; Proctor et al. 
2011). The high NFW (Navarro et al. 1997) halo con- 
centration parameters, lack of spectral star formation 
indicators, and large AR magnitude difference suggest 
these systems finished merging in the distant past, per- 
haps before z ~ 1 (Jones et al. 2000; Wechsler et al. 2002; 
D'Onghia et al. 2005; Khosroshahi et al. 2007). 

Recent studies of X-ray selected FGs paint a more com- 
plicated picture. The cooling time of FGs is significantly 
shorter than the Hubble time (e.g., Sun et al. 2004; Khos- 
roshahi et al. 2004, 2006a), yet they typically lack cool 
cores, suggesting that these systems may be younger or 
more active than previously thought (Mendes de Oliveira 
et al. 2009). Regular (non- fossil) rich groups often pos- 
sess cool cores (e.g., Finoguenov & Ponman 1999), even 
in the presence of AGN activity. In addition, there is 
evidence for enhanced SN II metal fraction in the cen- 
tral regions of FGs, suggesting a scenario where SN II 
powered winds resulting from merging late type galaxies 
erase the original central SN la Fe mass fraction domi- 
nance (Dupke et al. 2010). This is consistent with the 
previously found disky isophotes of the central dominant 
galaxies in FGs by Khosroshahi et al. (2006b) and also 
with the presence of shells in the stellar component in 
at least one of these galaxies, indicative of multiple past 
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mergers (Eigenthaler & Zeilinger 2009). Furthermore, 
some authors suggest that the FGs we see are the tail 
of the cluster distribution, possessing few galaxies at 
their current epoch for any number of reasons: failure to 
form those galaxies, early merging, or a quiescent state 
during a cycle of galaxy accretion (Mulchaey & Zablud- 
off 1999; von Benda-Beckmann et al. 2008; La Barbera 
et al. 2009; Dariush et al. 2010; Smith et al. 2010; Cui 
et al. 2011). The truth about fossil groups is somewhat 
muddled by the phenomenological rather than physical 
definition of the class. Extracting a useful physical defini- 
tion is in turn complicated by the relatively small number 
of FGs with deep X-ray observations. Indeed, what we 
call "fossil groups" perhaps comprise a heterogeneous set 
of galaxy systems with different formation and evolution 
histories. 

To characterize the ages and structural properties of 
FGs, it is crucial to have good data, especially X-ray 
observations. Detailed study of the intra-group medium 
(IGM) metal abundance, temperature structure, and in- 
ferred mass distribution help to constrain the halo forma- 
tion epoch and the importance of recent star formation 
or AGN activity. The available X-ray data are typically 
photon-poor due to the serendipitous nature of FG de- 
tections, and this has limited their study. To address 
this problem, we have embarked on a project to identify 
a large sample of FGs for future detailed follow-up stud- 
ies. In the work presented here, we have constructed a 
sample of 15 fossil group candidates, using the maxBCG 
cluster catalog (Koester et al. 2007a) to optically identify 
the candidates, which are then targeted with X-ray snap- 
shots using the Chandra X-ray Observatory to confirm 
the existence of a bright X-ray halo. This initial sample 
of confirmed FGs (using the phenomenological classifi- 
cation) will be invaluable for follow-up, including deep 
X-ray observations to study the metallicity structure of 
the gas and concentration of the mass distribution; and 
optical spectroscopy to compare the velocity dispersion 
to the X-ray mass. 

Throughout this paper we use a ACDM cosmology 
with ft m 0.27, tt A 0.73, and H 70 km s" 1 Mpc" 1 
(or h = 0.7). Except where specified with h notation, all 
numerical values from the literature have been scaled to 
correspond to this cosmology. Uncertainties are la and 
upper/lower limits are 3<r, unless stated otherwise. 

2. OPTICAL SAMPLE SELECTION 

To select a large sample of FG targets, we used the 
maxBCG cluster catalog (Koester et al. 2007a). This is 
a volume-limited catalog of over 17,000 optically selected 
red-sequence clusters in the redshift range 0.1 < z < 0.3 
with precise photometric redshift s (5 Z ~ 0.01) and op- 
tical richness estimates from the Sloan Digital Sky Sur- 
vey DR4 (SDSS; York et al. 2000; Adelman-McCarthy 
et al. 2006). The optical richness employed, A^oo, is the 
number of red-sequence member galaxies brighter than 
0.4 L* (in the z-band) found within a scale radius r2oo of 
the brightest cluster galaxy (BCG; Hansen et al. 2005). 
The mean properties of this catalog have been studied 
in detail, and we have obtained mean X-ray luminosi- 
ties (Rykoff et al. 2008), velocity dispersions (Becker 
et al. 2007), and a mass calibration via the mean weak- 
lensing shear profile around maxBCG clusters (Johnston 
et al. 2007; Sheldon et al. 2009). Simulations have shown 



that the catalog purity and completeness are very high 
(>90%; Koester et al. 2007b; Rozo et al. 2007). 

We selected for several optical characteristics that are 
expected of FGs from the empirical definition of the class. 
In particular, at a given optical richness, these systems 
should have a larger magnitude difference between the 
BCG and the next brightest galaxy and should have 
highly luminous BCGs typical of massive clusters. In 
our initial selection we restricted ourselves to systems in 
the richness range 9 < A200 < 25, corresponding to a 
mass range of 3xl0 13 < M 200 < lxlO 14 /*" 1 M© (John- 
ston et al. 2007). Our aim was to select systems rich 
enough to have sufficient X-ray luminosity for detection 
and analysis, and at the same time remain in the FG 
range where the BCG can truly dominate the system. 
We note that the publicly released maxBCG catalog was 
restricted to A^oo > 10; we used a slightly extended cat- 
alog described by Rykoff et al. (2008), including an addi- 
tional 3532 clusters with 7V 2 oo = 9 and allowing us to use 
the same richness bins studied in the analysis of maxBCG 
galaxy dynamics (Becker et al. 2007) and weak lensing 
(Johnston et al. 2007; Sheldon et al. 2009). To maximize 
the flux in the X-ray, we restricted the sample to those 
systems with a confirmed spectroscopic BCG redshift in 
the range 0.09 < z < 0.15; the lower redshift cut is im- 
posed by the maxBCG photometric redshift lower limit 
of z = 0.10 ± 0.01. At the time of selection from SDSS 
DR4, 42% of all maxBCG clusters in this redshift range 
had a spectroscopically determined BCG redshift. 

To quantify the galaxy magnitude gap representative 
of FGs, we used the difference in i-band magnitude of 
the BCG and the next brightest red-sequence cluster 
member within 0.5r2oo> denoted by A$. The r2oo values 
were estimated from the mass-scaling relation of John- 
ston et al. (2007). Bright non-red-sequence galaxies pro- 
jected within 0.5r2oo were not considered because the 
majority of these are foreground galaxies unassociated 
with the cluster. Thus, we have opted for a more com- 
plete sample of systems with this magnitude gap at the 
risk of a small amount of impurity. We refined our selec- 
tion to all the systems with A^ > 2.0 and BCG i-band lu- 
minosity Lbcg > 9xl0 10 /*- 2 L© (Lbcg > 1.8xlO n L© 
in our adopted cosmological framework). There were 26 
maxBCG systems that passed this selection cut. We il- 
lustrate the maxBCG selection in Figure 1, plotting A^ 
as a function of Lbcg] these two values are correlated, as 
the most luminous BCGs also tend to have a large A^. 
However, richer maxBCG clusters (A/200 > 25) in the 
same redshift range (magenta points) do not extend this 
trend to the largest A;, despite having relatively larger 
Lbcg- We discuss selection effects of our sample more 
fully in a companion paper by Proctor et al. (2011). The 
cuts based on these parameters are consistent with the 
majority of known FG systems (e.g. Khosroshahi et al. 
2007) that overlap the SDSS footprint, marked with blue 
diamonds in Figure 1. These five systems are present in 
the maxBCG catalog, however they each fail one or more 
of our selection criteria, with four falling out of our red- 
shift range and the fifth having A/200 = 28. 

Deeper X-ray follow-up is a major goal of this sam- 
ple assembly, and the most efficient current instrument 
for this purpose is the XMM-Newton X-ray Observa- 
tory. The ~ 15 arcsec spatial resolution of XMM-Newton 
corresponds to 30 kpc at z = 0.1. This is a sizeable 
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Fig. 1. — Difference in i-band magnitude (A^) of the BCG and 
next brightest red-sequence cluster member within 0.5r2oo as a 
function of BCG i-band luminosity (Lbcg) f° r an maxBCG sys- 
tems with measured spectroscopic BCG redshifts 0.09 < z < 0.15. 
Our FG candidates (red squares) were chosen from the subset of 
systems with 9 < ./V200 < 25 (black dots) within the region de- 
lineated by the dotted lines; for comparison we also show all sys- 
tems with A^2oo > 25 (open magenta circles). Known FGs from 
Khosroshahi et al. (2007) that overlap the SDSS footprint and are 
detected by maxBCG (blue diamonds) are shown for illustration 
although they do not meet our selection criteria. Note that three 
of our candidate FGs have very similar values for A* (2.25-2.26) 
and Lbcg (10xl0 10 L©), so the red squares overlap. 

fraction of the typical group core radius, so to maxi- 
mize the utility of these targets for deeper X-ray spec- 
troscopic analysis we sought to reduce the possibility 
of bright central AGN emission that might contaminate 
the diffuse IGM emission. We rejected all BCGs that 
have evidence of Seyfert or LINER-like line emission with 
log([NII]/Ha) > -0.2 (e.g. Kauffmann et al. 2003). This 
eliminated 3 of the 26 systems. Additionally, we rejected 
all BCGs that match radio sources in the FIRST catalog 
(White et al. 1997) within 3 arcsec of the BCG position 
(6 kpc at z = 0.1), eliminating another 6 candidates. Fi- 
nally, we rejected 2 systems that have bright stars with 
ttir < 7 in the field-of-view that exceed the bright star 
limit for XMM-Newton. In total, we identified 15 can- 
didate FGs, marked with red squares in Figure 1 and 
summarized in Table 1. 

The elimination of candidates based on AGN activ- 
ity, especially in the radio, could produce a selection 
bias. The motivation for this selection was based on 
Allen et al. (2006) and Balmaverde et al. (2008), who 
find a clear correlation between Bondi accretion rate and 
central engine jet power in samples of low-power radio 
galaxies. The presence of radio emission from a jet or 
expanding bubbles is a signal of strong accretion, with 
rate MBondi c£ P- As the X-ray emission measure of this 
plasma is EM oc p 2 , the flux within the central few kpc 
can be quite bright and contaminate the XMM-Newton 
surface brightness profile, and through an abundance of 
caution, we expunged these radio-bright candidates. The 
direct correlation between AGN radio and X-ray lumi- 
nosity in BCGs is not well constrained, although Hickox 
et al. (2009) find a small overlap (<~ 10%) between radio- 
and X-ray-bright AGN in general, with the former tend- 
ing to live in luminous red sequence galaxies and the 
latter tending to "green valley" galaxies. The possible 
effects of this selection on our sample are discussed in 
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Fig. 2. Radial stacked RASS profile of 170 maxBCG clusters 
with a similar selection function to those targeted in this paper. 
The thick dashed line shows the best-fit /3-model surface brightness 
profile, which yields a good fit with (3 = 0.6 ±0.1. The dotted lines 
show profiles for a RASS point source with two normalizations; 
the lower one is scaled to the central surface brightness, while the 
higher one is scaled to the surface brightness at 0.4 h -1 Mpc. In 
both cases it is clear that the stacked X-ray emission of this sample 
is significantly extended when compared to a RASS point source. 



Section 4.1. 

The final criterion in the FG definition is a cluster- 
scale X-ray halo, which also serves as confirmation of a 
collapsed system rather than a projection of unrelated 
galaxies. While six of the FG candidates are bright 
enough to be detected at the > 2a level in the ROSAT 
All Sky Survey (RASS), a majority of these optically se- 
lected targets do not have sufficient flux to be detected in 
RASS (see Table 1). We therefore estimated the mean 
Lx and radial profile using the stacking procedure de- 
scribed by Rykoff et al. (2008), selecting 170 maxBCG 
systems chosen with criteria used for the 15 FG targets 
(9 < A^200 < 25; A* > 2.0; and L B cg > 9 x 10 10 L©). In 
order to select enough systems for the stacking analysis 
we expanded the redshift range to 0.09 < z < 0.20, and 
we did not require BCG spectroscopic redshifts nor did 
we filter for active galaxies or bright stars. The stacked 
RASS profile is shown in Figure 2 along with a f3 model 
fit with p = 0.6 ± 0.1. To show that the stacked FG 
profile is significantly extended, we also stacked a rep- 
resentative sample of RASS point sources treated as if 
these sources were at the redshifts of the maxBCG tar- 
gets (see Section 3.3.1 in Rykoff et al. 2008). The stacked 
point source profile is shown with dotted lines with two 
scalings in Figure 2. In one scaling the central flux is 
matched to the stacked FG profile; in the second, the 
flux at 0.4 Mpc is matched to the FG profile. In both 
cases it is clear that the stacked FG emission is signifi- 
cantly extended when compared to RASS point sources. 

The derived Lx estimates were used to plan our follow- 
up X-ray observations with Chandra, described in the 
next Section. Rykoff et al. (2008) show that the mean 
Lx of maxBCG systems scales as a power-law with A^oo 
(over two orders of magnitude in Lx) and also scales 
with Lbcg fc> r the poorer clusters and groups. Thus, 
the FGs are expected to be more X-ray luminous than 
typical maxBCG systems at a similar richness. 
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TABLE 1 

Fossil Group Sample and Chandra Observations 





Target BCG 




short name a 


RA b 


Dec b 




z c 


A^200 d 


(mag) 


Lbcg { 
(10 10 L©) 


RASS§ 


obs. date (OBSID) 


/ h 

^exp 

(ksec) 


SDSS 


J013325 


87-102618 


6 


J0133-1026 


23 


3578 


-10 


4385 





113 


12 


2.45 


10 


4 




2009-05-29 


(10753) 


10.0 


SDSS 


J081526 


59+395935 


5 


J0815+3959 


123 


8608 


+39 


9932 





129 


12 


3.26 


16 


7 


V 


2008-12-12 


(10758) 


5.1 


SDSS 


J082122 


54+405123 


7 


J0821+4051 


125 


3439 


+40 


8566 





125 


10 


2.19 


9 


9 


2009-01-04 


(10474) 


10.0 


SDSS 


J085640 


72+055347 


3 


J0856+0553 


134 


1697 


+5 


8965 





094 


16 


2.26 


10 







2009-01-09 


(10750) 


5.5 


SDSS 


J090638 


27+030139 


1 


J0906+0301 


136 


6595 


+3 


0276 





136 


9 


2.93 


10 


9 




2009-01-14 


(10475) 


10.0 


SDSS 


J100742 


53+380046 


6 


J1007+3800 


151 


9272 


+38 


0130 





112 


24 


2.54 


16 


9 


V 


2009-02-09 


(10755) 


4.7 


SDSS 


J101745 


57+015645 


8 


J1017+0156 


154 


4399 


+1 


9461 





118 


12 


2.34 


11 


2 




2009-03-23 


(10754) 


9.9 


SDSS 


J103930 


43+394718 


9 


J1039+3947 


159 


8768 


+39 


7886 





093 


14 


2.46 


9 


9 




2009-01-14 


(10749) 


5.1 


SDSS 


J104548 


50+042032 


5 


J1045+0420 


161 


4521 


+4 


3424 





154 


13 


2.07 


9 


6 




2009-02-01 


(10476) 


9.9 


SDSS 


J113305 


51+592013 


7 


J1133+5920 


173 


2730 


+59 


3372 





133 


13 


2.26 


10 





V 


2009-07-08 


(10472) 


5.7 


SDSS 


J113623 


71+071337 


5 


J1136+0713 


174 


0988 


+7 


2271 





103 


17 


2.25 


10 





V 


2009-02-09 


(10756) 


5.0 


SDSS 


J115305 


32+675351 


5 


J1153+6753 


178 


2722 


+67 


8977 





117 


17 


2.19 


13 


4 


V 


2009-06-21 


(10473) 


5.0 


SDSS 


J133626 


96+545353 


8 


J1336+5453 


204 


1124 


+54 


8983 





107 


10 


3.02 


13 


9 


2009-09-25 


(10752) 


7.1 


SDSS 


J141004 


19+414520 


8 


J1410+4145 


212 


5175 


+41 


7558 





094 


21 


2.30 


12 


7 


V 


2009-07-07 


(10757) 


5.1 


SDSS 


J141115 


89+573609 





J1411+5736 


212 


8162 


+57 


6025 





106 


16 


2.19 


9 


2 


2009-07-23 


(10751) 


6.9 



a The short names for each target are used throughout this work. 
b RA, Dec are the J2000 coordinates of the BCG, in degrees. 

c The redshift is the BCG spectroscopic value from SDSS DR4 (Adelman-McCarthy et al. 2006). 

d A^200, a richness estimate, is the number of red-sequence cluster galaxies brighter than 0.4 L* (in the i-band) found within r2oo of the BCG. 
e Ai is the difference in i-band magnitude of the BCG and the next brightest red-sequence cluster member within 0.5r200- 
f Lbcg is the i-band luminosity of the BCG. 

g RASS indicates a 2a or better detection in the ROSAT All Sky Survey. 
h Effective exposure time of cleaned event data. 



3. X-RAY OBSERVATIONS AND DATA ANALYSIS 

3.1. X-Ray Data and Reduction 

The X-ray observations were performed with the Chan- 
dra X-ray Observatory between December 2008 and 
September 2009 in the form of 5-10 ksec snapshots (see 
Table 1). Data were obtained using the ACIS-S3 chip, 
with the candidate fossil group BCG centered in the field 
of view. Standard processing was performed on the raw 
event files, including the background reduction tools ap- 
plicable to the VFAINT observing mode 13 . The resulting 
0.3-7 keV light curves were filtered to remove additional 
times of high background and applied to produce cleaned 
event files. Point sources were identified from the 0.3-7 
keV events using the CI AO tool wavdetect and masked 
out for the analysis of the extended emission. This in- 
cluded any emission clearly identified with the optical ex- 
tent of the BCG, whether point-like or slightly extended. 
A separate analysis of these features is presented in Sec- 
tion 3.4. 

3.2. X-Ray Spectral Analysis 

Spectral extraction regions were identified from images 
in the 0.5-2 keV band, where the group emission should 
dominate, and chosen to encompass the bulk of the ex- 
tended X-ray emission, with radii in the range X.I 1 — 2. 1 1 
(129-356 kpc) centered on the peak of the X-ray emis- 
sion. For two targets with irregular X-ray morphology 
(J0133-1026 and J1045+0420), the region was centered 
on the apparent centroid of the emission. Targets with- 
out obvious group emission (J0821+4051, J0906+0301, 
and J1336+5453) were assigned an extraction region of 
250 kpc (1.7 7 — 2.1 7 ) centered on the optical BCG location 
to estimate an upper limit on the flux. A background 

13 See http://cxc.harvard.edu/cal/Acis/Cal_prods/bkgrnd/ 
current /background . html. 



region was defined for each observation from the remain- 
ing area of the 8.4' x 8.4' ACIS-S3 chip, excluding the 
inner ~ 2.5' radius, the identified point sources, and the 
outer edge of the field of view. The outer extent of each 
background region was typically 7.5' x 5.5' in size, ori- 
ented along the direction of the ACIS-S CCD array. This 
choice was made to reduce the effects of non-uniform 
molecular contamination on the CCD, which is thought 
to be thicker near the edges of the ACIS-S array and 
which substantially reduces the soft X-ray transmittance 
(Vikhlinin 2004). 

The spectral analysis was performed in XSPEC vl2.6.0 
utilizing the C-statistic, a modified Cash (1979) likeli- 
hood function that allows for inclusion of a background 
spectrum and a goodness-of-fit estimator similar to x 2 
in the limit of many counts 14 . The X-ray spectrum for 
each group (see Figure 3) was fit with an absorbed APEC 
model (Smith et al. 2001) in the 0.4-7 keV band, with the 
redshift fixed at the BCG spectroscopic value from SDSS. 
The intervening Galactic TVhi column was fixed at the 
average value reported by the Leiden/ Argentine/Bonn 
(LAB) merged survey (Kalberla et al. 2005); in all cases 
it was less than 5 x 10 20 cm -2 . The temperature and 
normalization were allowed to vary. For several groups, 
the metal abundance (using Anders & Grevesse (1989) 
photospheric solar abundances) was unconstrained dur- 
ing the initial fit and was frozen at 0.3 solar, the weighted 
mean value from the well-constrained fits and similar to 
that measured in ~ 2 keV systems (Osmond & Ponman 
2004; Rasmussen & Ponman 2007). The spectral fitting 
results are shown in Table 2. 

We consider a detection to be a 3a or greater excess of 
counts in the 0.5-2 keV band compared to the expected 

14 The XSPEC implementation of the C-statistic is described 
in detail at http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 
manual/XSappendixCash . html. 
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Fig. 3. — Background-subtracted Chandra/ ACIS-S3 spectra of the 12 detected FGs with the best-fit models. The displayed spectra have 
been binned in energy for clarity, although the spectral analysis was carried out on the full resolution (14.6 eV/channel) ACIS spectra. 
Shown below each spectrum are the fit residuals. The FGs are shown roughly in order of brightest to faintest for clarity of the ordinate 
scaling. 



background within the extraction region. Twelve of the 
15 targets were detected using this definition; all of the 
detections are easily visible to the eye in point-source- 
excluded images (see Figure 4). Eleven of these detected 
systems have well-constrained temperatures in the range 
1-3 keV; while J0133-1026 is detected at 4.9a, the 104 
source counts are insufficient to constrain the spectral 
model. Errors on the spectral parameters were deter- 
mined by sampling parameter space for each parameter, 
marginalizing over the other free parameters. For the 
systems with fixed abundance, we estimated the errors 



in other parameters by stepping (with XSPEC steppar) 
the abundance over the expected range of 0.1-1 solar. 
Absorbed fluxes and unabsorbed "soft" (0.5-2.0 keV rest 
frame) luminosities were determined from the best-fit 
spectral models. A "bolometric" (0.008-100 keV rest 
frame) luminosity with errors was determined for each 
group by extrapolating the unabsorbed spectral model. 
We assumed kT = 2 keV and abundance of 0.3 solar for 
the three undetected targets, and used this model along 
with the background counts estimate to determine up- 
per limits on the diffuse X-ray flux and luminosity. An 
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Fig. 4. — SDSS multi-color images of the 12 CTiandra-detected FGs, with the diffuse 0.5-2 keV X-ray emission overlaid in contours. The 
Chandra contours have point sources removed and are exposure-corrected and smoothed with a 30" FWHM Gaussian. The displayed blue 
intensity and contours start at 2a (about 5xl0 -9 ph s — 1 cm -2 ) above the background surface brightness, and increase in intervals of 2a. 
Red circles note the optically identified BCG, and the green circle indicates 0.5r5oo, centered on the X-ray peak (for well-defined /3-model 
fits) or BCG position. The green bar in the upper left of each panel shows 100 kpc at the FG redshift. Images are 6' on a side, ranging 
from 0.62 to 0.97 Mpc for the nearest and farthest FG, respectively. 



identical model was assumed for J0133— 1026 to calcu- 
late its flux and luminosity. The source and total counts 
are shown in Table 2, along with the flux and luminosity 
estimates. 

3.3. X-Ray Spatial Analysis 

To compare to existing results and expectations for 
self-similar scaling, cluster X-ray luminosities and masses 
are typically scaled to a common radius in terms of the 
average interior overdensity S r = p r /pcrit, where p r is 
the mean cluster mass density within radius r, and p cr i t 
is the critical density at z (e.g., Maughan et al. 2006). 
The extraction regions we have applied are considerably 
smaller than the typical radius for S r = 500 for group 
potentials (Y500 ^ 500 kpc, compared to r extract ~ 250 
kpc; see Table 3). Therefore we must estimate the lumi- 
nosity correction factor by extrapolating the X-ray sur- 
face brightness profile, similar to the approach in previ- 
ous studies (e.g., Jeltema et al. 2006). In the following 
analysis, we use the working assumption that the intra- 
group medium is spherical, non-rotating, isothermal, and 
in hydrostatic equilibrium to r$oo within the group grav- 
itational potential. 



Counts images of the ACIS-S3 field of view were con- 
structed from the cleaned event lists, binning to 4 x 4 
pixels (2x2 arcsec) and restricting the energy band to 
0.5-2 keV, where the group emission dominates the back- 
ground. For each detected group, two-dimensional spa- 
tial fitting was performed with the Sherpa package avail- 
able in CI AO, using a circular f3 model surface brightness 
profile for the FG emission and a constant baseline to 
account for the combined cosmic and instrumental back- 
ground. The j3 model was multiplied by an exposure 
map during the fit, while the background was multiplied 
by a mask containing the bad pixels and columns of the 
CCD, dithered according to the aspect solution of the 
observation. Fits were performed in two dimensions to 
a region ~ 7.5' x 5.5' oriented along the direction of the 
ACIS-S array, encompassing the region defined as the 
background for spectral fitting (see Section 3.2), chosen 
in such a way as to to reduce the effects of non-uniform 
molecular contamination on the CCD. Point source re- 
gions were excluded in the fitting, which was done using 
the Cash (1979) statistic, allowing /3, core radius r c , the 
emission center, and the FG and background amplitude 
to vary. The fit results are summarized in Table 3. 
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Spectral Fitting Results 
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a kT was fixed to 2 keV to estimate detection limits for undetected sources and J0133— 1026. 

b Abundance assumes the solar photospheric values of Anders Sz Grevesse (1989). The value was fixed to 0.3 for fits 
with unconstrained abundance. 

c Absorbed model flux in the 0.5-2 keV band. Upper limits are 3cr. 

d Unabsorbed model luminosity in the 0.5-2 keV band, rest frame. Upper limits are 3cr. 

e Unabsorbed model luminosity in the 0.008-100 keV band, rest frame. Upper limits are 3cr. 

f Counts in the spectral extraction region in the 0.5-2 keV band, observed frame. Source counts are estimated from 

the spectral model. 

g Detection significance in units of a. 
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Spatial Fitting Results 
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a Aperture correction factor to convert Lj through the observed aperture to rsoo- 

b Unabsorbed model luminosity in the 0.5-2 keV band, rest frame, corrected to rsoo- 

c Unabsorbed model luminosity in the 0.008-100 keV band, rest frame, corrected to rsoo- 
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Of the 12 FGs detected in extended X-ray emission, 
four have sufficient counts to constrain the f3 model pa- 
rameters. These are four of the five brightest targets, 
with F > 20xl0- 14 erg cm" 2 s" 1 (0.5-2 keV). The fifth 
bright target, J1045+0420, is morphologically irregular 
and obviously not well-fit by a simple f3 profile (see Sec- 
tion 4.3). The best-fit f3 values range from 0.3 to 0.5, 
smaller than the value of 0.67 commonly found for clus- 
ters (e.g., Jones & Forman 1999) but not unusually small 
for rich groups and poor clusters of similar temperature 
in this redshift range (Willis et al. 2005; Jeltema et al. 
2006) or at z ~ (Osmond & Ponman 2004). The values 
for the core radius r c are all comparatively small, ranging 
from 7 to 50 kpc, but consistent with the previously cited 
results. For the seven detections with unconstrained f3 
model parameters, we assumed the average best-fit val- 
ues of fi = 0.4 and r c = 25 kpc to estimate the spatial 
extent of the X-ray emission and the luminosity correc- 
tions. Note that despite an X-ray detection, J0133— 1026 
was excluded from the remaining analysis, since its lack 
of a measured temperature rendered the spatial extrap- 
olation too uncertain. 

Based on the best-fit /cT, /?, and r c with associated 
errors, we estimated rsoo and M 5 qo for each group. The 
mass within radius r can be given as (e.g., Arnaud & 
Evrard 1999): 

M(< r) = 1.13X10 14 P^-Z- - {V { r j\ 9 M . (1) 
v ) H keV Mpc 1 + (r/r c ) 2 u v ' 

Since rsoo ^> r ci we can estimate rsoo with a simple 
analytic approximation (e.g., Arnaud & Evrard 1999, 
Eq. 17). We chose to iteratively solve Eq. 1 for M500 
and r5oo, using the definition 



TABLE 4 
BCG X-ray Emission 
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The results are shown in Table 3; r^oo = 443-712 kpc 
and M500 = 0.3-1.0xl0 14 M for the sample, typical 
values for groups and clusters in this temperature and 
redshift range (Willis et al. 2005; Jeltema et al. 2006; 
Finoguenov et al. 2007; Jeltema et al. 2009). These val- 
ues are equivalent to what we derive from the analytic 
approximation. 

Aperture corrections were calculated to scale the ob- 
served luminosity to that within r$oo , including the small 
(< 2%) correction for excluded point source regions. The 
aperture corrections range from 1.5 to 3.0 with an aver- 
age correction of 2.2. While these are large corrections, 
the attendant errors take into account the uncertainty 
in the j3 model parameters. The corrections for the four 
groups with well-constrained f3 models are 1.6-2.8, with 
errors of less than 30%. The aperture corrections and 
corrected r$oo luminosities are listed in Table 3. 

To analyze the structure of the hot intra-cluster 
medium (ICM), researchers typically employ power ra- 
tios (multipole moments of the X-ray surface brightness) 
and centroid shifts that are sensitive to substructure and 
irregular morphology (e.g., Jeltema et al. 2008). Power 
ratios require thousands of counts and a well-constrained 
surface brightness out to ~ r^oo , therefore we cannot em- 
ploy them with the current snapshot data. We attempted 
to calculate centroid shifts for the detected FGs, however 
the results were inconclusive, with large errors driven 
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a T is the power law index, fixed to 1.7 for all but one source. 
b Absorbed model flux in the 0.5-2 keV band. 
c Ratio of the BCG 0.5-2 keV flux to the detected diffuse flux from 
Table 2. 

d Unabsorbed model luminosity in the 0.5-2 keV band, rest frame. 
e Counts in BCG region, 0.3-7 keV. No background has been 
subtracted. 

by the low counting statistics. We present a qualitative 
analysis of the morphology, specifically the fraction of 
disturbed and relaxed clusters, in Section 4.1. 

3.4. BCG X-ray Emission 

To reduce the possibility of contamination from AGN 
or diffuse emission in the BCG, we excluded sources 
found by wavdetect in Section 3.1 that fall within the 
r-band optical extent of the BCG. As the images in Fig- 
ure 5 show, 8 of the 12 BCGs possess a detected central 
source with emission in the 0.3-7 keV band, ranging from 
7 to 35 counts in the regions outlined by red ellipses. Sev- 
eral of these sources appear to be point-like AGN, but 
all 12 FGs have some extended component as well, either 
from the BCG, the core of the IGM, or some combination 
thereof. 

We estimated the contribution of these detected 
sources to the total X-ray flux by extracting a spectrum 
for each source and fitting with XSPEC. We assumed a 
spectral model typical of an AGN, an absorbed power law 
with fixed power law index T = 1.7, and used the fixed 
Nh and redshift values from Table 2. The power law nor- 
malization was allowed to vary, which provided a confi- 
dence interval for the source flux. For one of the targets, 
the source in J0133— 1026 with 35 counts, we allowed T 
to vary and obtained a good constraint of T = 1.8^0.4- 
The 0.5-2 keV flux estimates (see Table 4) range from 
a few percent to 30 percent of the flux from the diffuse 
IGM emission. 

One or more of the detected central X-ray sources 
could be thermal in origin, along the lines of the em- 
bedded coronae identified in cluster galaxies by Sun et al. 
(2007). This is especially true of the more extended emis- 
sion sources, such as those in SDSS J1007+3800, SDSS 
J1410+4145, and SDSS J1411+5736. While the spectral 
model used here is inappropriate for thermal BCG X-ray 
emission, it nevertheless provides a reasonable estimate 
of the flux in the observed band for planning future deep 
X-ray observations. Due to the paucity of counts from 
these sources, a detailed spectral analysis is not possi- 
ble with the current data. FGs with notable BCG X-ray 
flux are discussed in more detail in Section 4.3. We em- 
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Fig. 5. — Chandra 0.3-7 keV counts image (left) and SDSS r-band image (right) for each FG, centered on the BCG. Blue dashed contours 
show the same point-source-excluded, smoothed emission as in Figure 4. Red ellipses identify BCG emission that has been excluded from 
the spectral and spatial analysis of the IGM emission for eight of the FGs, and analyzed separately in Section 3.4. Images are 30" on a 
side. 



phasize that we have excluded the BCG emission in our 
spectral (Section 3.2) and spatial (Section 3.3) analysis 
for the 8 clusters in which we detect it. Therefore the 
total X-ray flux for an instrument which is unable to re- 
solve and exclude these central sources would be the sum 
of the values listed in Tables 2 and 4. 

4. RESULTS AND DISCUSSION 

4.1. Success of the Method: Are They 
"Real" Fossil Groups? 

The primary goal of this work, to optically identify a 
sizable sample of fossil groups, has met with great suc- 
cess. Out of a sample of 15 candidates, we have con- 
firmed 12 gravitationally bound systems through their 
IGM X-ray emission. One additional target, SDSS 
J0906+0301, is undetected in the short Chandra expo- 
sure but is clearly a bound system from follow-up optical 
spectroscopy, with a velocity dispersion of a = 506 ± 72 
km s _1 based on 25 member redshifts (Proctor et al. 
2011). This 87% success rate demonstrates the value of 
the maxBCG survey for selecting FGs; moreover, since 
these galaxy-poor systems are the most difficult to iden- 
tify optically, we expect a very high success rate for this 
method in selecting clusters and groups in general. How- 
ever, we note two caveats that must be addressed. 

First, our original FG criteria were based on a mass- 
scaling relation from SDSS systems (Johnston et al. 
2007). From the current analysis and the results of Proc- 



tor et al. (2011), the systems in our sample are very 
massive for their richness, and in fact they are under- 
luminous in the optical (including the BCG luminosity) 
by about a factor of three. Simply scaling for the ob- 
served richness could dramatically underestimate r2oo? 
which in turn could exclude bright galaxies from our A$ 
magnitude difference criterion and render our systems 
non- fossil groups. To address this, we have calculated 
new values for r2oo based on the X-ray results presented 
here. In addition to the the f3 model approach with which 
we estimated rsoo, we adopt the scaling relation of Hels- 
don & Ponman (2003), 



^200, kT 



1.14 



kT 
keV 



1/2 



E(z)- 1 Mpc, (3) 



where E(z) = H(z)/H = [ft m (1 + zf + Ov] 1 / 2 for 
a ACDM universe (e.g., Maughan et al. 2006). These 
values and the j3 model extrapolation values (denoted 
^200, p) are presented in Table 5, along with the magni- 
tude difference determined for each r2oo and its 1<j 
spread of values. The r2oo values originally used for the 
maxBCG optical selection are hereafter denoted r2oo,m#- 
The original r2oo,m# values are consistent within the 
errors of the f3 model results for all but two FGs. For 
J1410+4145, r2oo,m£ was an overestimate, so no addi- 
tional bright galaxies are included using the r2oo,/3 value. 
For J 10 17+0 156, r2oo mB was an underestimate, and one 
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TABLE 5 

Comparison of r2oo Estimates and Resulting A* 
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(2.3,0.5) 


1266 (1123,1608) 


0.5 (0.5,0.5) 


J1153+6753 


17 


855 


2.2 


825 (661,1041) 


2.2 


(2.2,2.2) 


1023 (892,1295) 


2.2 (2.2,1.2) 


J1410+4145 


21 


934 


2.3 


784 (720,869) 


2.3 


(2.3,2.3) 


997 (937,1084) 


2.3 (2.3,2.3) 


J1411+5736 


16 


833 


2.2 


783 (584,978) 


2.2 


(2.2,2.2) 


974 (797,1221) 


2.2 (2.2,1.6) 



a A^2oo is defined in Table 1. 

b ^200 estimate from Johnston et al. (2007) used in our maxBCG optical selection, as discussed in Section 2. 
c A^ is defined in Table 1. Each column is determined from the r2oo preceding it, with values in parentheses 
from the la extrema of r2oo- 

d T200 estimate from extrapolating our (3 model fit, with la extrema in parentheses. 

e ^200 estimate from Helsdon &; Ponman (2003) (also see Eq. 3), with la extrema in parentheses. 



bright galaxy falls within r2oo,/3 such that A$ = 1.8 mag. 
This barely fails the magnitude difference criterion of 
A^ > 2 mag; in the r-band, the difference is only 1.9 
mag (Proctor et al. 2011). Two other FGs (J1039+3947 
and Jl 136+0713) are consistent with the A^ criterion at 
the lower bound of r2oo,/3, but inconsistent at its upper 
bound. 

The r2oo,fcT values are systematically ~ 20% higher 
than r2oo,mB and r2oo,/3, however, the confidence inter- 
vals overlap for 8 out of the 11 FGs. Only the uncertainty 
in the X-ray kT measurement has been propagated to 
the r 200, kT confidence estimates, therefore any intrinsic 
scatter in the Eq. 3 scaling relation would increase the 
uncertainty. With the current confidence intervals, two 
FGs (J1017+0156 and J1136+0713) fail the A* criterion 
for all ±la values of r 200, kT- These systems are noted in 
Figure 6 with blue squares. A total of 8 FGs fail that cri- 
terion for the extreme +la value, but nevertheless agree 
with the criterion within the uncertainty. Proctor et al. 
(2011) reach identical conclusions for the subset of 5 of 
these FGs for which they have r-band photometry. 

The second caveat concerns our X-ray surface bright- 
ness extrapolation. The j3 and core radii we have used 
are somewhat smaller than what is typically found for 
clusters (f3 ~ 0.67, r c ~ 100-200 kpc; Jones & Forman 
1999). Indeed, our extraction radii are not much larger 
than the typical r c value, and it is possible the low f3 
values result from fitting a fairly fiat core of emission. 
For j3 = 0.67 and our assumed r c = 25 kpc, the lumi- 
nosity aperture corrections to r$oo would be about 50% 
lower than what is listed in Table 3. For r c = 100 kpc, 
the aperture corrections are 25% lower, and for r c = 200 
kpc, they are 20% higher. These systematic uncertainties 
are within the statistical errors quoted for the aperture 
corrections. 

In summary, the A^ magnitude difference criterion is 
satisfied for the majority of our FGs with measured tem- 
perature, insofar as our knowledge of r2oo is correct, and 
therefore we conclude that these systems are real fossil 
groups. Deeper X-ray observations will more precisely 
constrain the kT and surface brightness profiles to al- 
low more accurate hydrostatic mass determinations. Re- 



gardless of the details of whether these systems strictly 
meet the empirical FG definition, the results presented 
by Proctor et al. (2011) clearly demonstrate that a sub- 
set of our FGs are different from other systems of similar 
mass. 

As noted in Section 2, the exclusion of candidates with 
known central radio sources could bias our sample. In 
particular, galaxy clusters with strong cool cores are very 
likely to have radio AGN (Sun et al. 2009; Mittal et al. 
2009). Thus to first appearances we have selected against 
strong cool core FGs. In fact, our screening criterion 
turns out to be insufficient to remove all radio sources, 
as 4 of the 12 confirmed FGs have radio lobes that are 
outside of the 3" radius cut that we imposed (see Sec- 
tion 4.3). Thus while some strong cool core FGs may 
have escaped our sample, it remains useful for future ra- 
dio studies, and we expect little bias in kT given the 
large measurement uncertainties in the current analy- 
sis. A deeper XMM- Newton study will allow removal 
of the core emission when assembling scaling relations, a 
technique that greatly reduces scatter and does not pro- 
duce any apparent bias (Pratt et al. 2009; Vikhlinin et al. 
2009a). 

Fossil groups typically have relaxed X-ray morphol- 
ogy, and the morphological differences in our sample 
are at first glance somewhat puzzling. As pointed out 
in Section 3.3, the modest signal-to-noise of our snap- 
shot data precludes a detailed analysis of the struc- 
ture of each FG. From a qualitative standpoint, we con- 
clude that much of the appearance of disturbed mor- 
phology is due to lack of photons; this is likely the 
case for J0133-1026, J0815+3959, J1039+3947, and 
J141 1+5736, and little can be said about the relaxed 
state of these four FGs. Five of the FGs are qualita- 
tively relaxed, having bright core regions and elliptical 
isophotes excluding the outer few contour levels (see Fig- 
ure 4); these are J0856+0553, J1017+0156, J1133+5920, 
J1153+6753, and J1410+4145. Of the three remaining 
FGs, J1045+0420 is clearly disturbed, J1007+3800 has 
irregular isophotes that are compressed on the east com- 
pared to the west, and Jl 136+0713 has broad but low- 
level irregular emission. Therefore, we can qualitatively 
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Fig. 6. — Lx~Tx relation for low- z groups and clusters, with the results overplotted for our 11 FGs with measured Tx- All Lx values 
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say that between 10% (1 out of 8) and 40% (3 out of 8) 
of the FGs are apparently not relaxed. 

These estimates for disturbed fraction are a bit lower 
than what is seen in observations and simulations of nor- 
mal systems. Bohringer et al. (2010) derive power ra- 
tios and centroid shifts for 31 clusters from the Repre- 
sentative XMM-Newton Cluster Structure Survey (REX- 
CESS) and identify two discrete samples in these param- 
eter distributions, with ~ 40% of the clusters identified 
as disturbed systems. This is in line with previous es- 
timates of the disturbed fraction (e.g., Jones & Forman 
1999). Jeltema et al. (2008) perform simulations of clus- 
ters and measure the same structure observables, and 
find that ~ 30% of the systems would appear disturbed 
in Chandra-quality X-ray images. Obviously the exact 
fraction of disturbed and relaxed clusters depends on 
how these classes are observationally defined and what 
metrics are used. From a qualitative view, the fraction 
of relaxed FGs in our sample is similar to what is seen 
and predicted in normal systems, and in fact could be 
significantly higher than that, more consistent with the 
expectation for FGs. 

4.2. Fossil Group Scaling Relations 

With the scaled luminosities, and with the caveats 
outlined in Section 4.1 in mind, we are able to place 
the FGs on a Lx~Tx relation. This is shown in Fig- 



ure 6 along with a number of low- z groups (Osmond & 
Ponman 2004) and clusters (Markevitch 1998; Wu et al. 
1999; Vikhlinin et al. 2009a; Pratt et al. 2009), with all 
values corrected with E(z) scaling for self-similar evolu- 
tion. All values have been scaled to h = 0.7 and are 
measured within or corrected to r$oo with the excep- 
tion of the Wu et al. (1999) clusters, which are measured 
within 1 Mpc. The FGs fall close to the locus of points 
traced by low- z systems, although they are consistent 
with rich groups and poor clusters instead of the poor 
group scales that one would infer from their richness. 
This is expected from previous observations of FGs (e.g. 
Khosroshahi et al. 2007) and indeed from the empirical 
definition of a fossil group. There is a hint from Fig- 
ure 6 that these FGs are either hotter or less luminous 
than normal systems, although the errors are large. The 
slightly larger aperture correction favored by a f3 = 0.67, 
r c = 200 kpc surface brightness model would move the 
FG points closer in Lx to the locus of points from normal 
systems. 

With our modest X-ray data we can explore other scal- 
ing relations that constrain different cluster characteris- 
tics. The Lx~M relation is a useful probe of total baryon 
fraction, which is dominated by the hot ICM. Using Eq. 1 
for r2oo, we calculate M200 for our FGs and extrapolate 
the X-ray luminosity to r2oo? using the method described 
in Section 3.3. The LX-M200 relation is shown in Fig- 
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Fig. 7. — LX-M200 relation from stacking of weak lensing measurements of maxBCG clusters (Rykoff et al. 2008b), with our FGs 
overplotted. The red star shows the geometric mean of our systems; it is clearly offset to lower mass and higher luminosity compared to 
the weak lensing results. However, systematic error in our /3-model surface brightness fitting is large enough to account for this offset, as 
demonstrated by the other star points, which show the centroid we would have obtained had we used values of j3 and r c typically found 
for normal clusters. 

ure 7, along with the relation from Rykoff et al. (2008b) 
for maxBCG clusters. We correct for self-similar evolu- 
tion at the redshift z = 0.25 used in this work, which uses 
stacked weak lensing measurements to estimate masses. 
Our FGs fall above the Rykoff et al. (2008b) fit, at higher 
luminosity for a given mass, which would indicate an en- 
hanced baryon fraction. However, the systematic error 
due to /3-model extrapolation (described above) is much 
larger than this typical offset; as illustrated in Figure 7, 
this systematic effect is more likely to underestimate the 
mass. The FGs are therefore consistent with Lx~M scal- 
ing relation of all maxBCG clusters, to the extent that 
we can constrain the hydrostatic masses of our FGs. 

Likewise, we compare the M-Tx relation to recent 
results from the literature (Vikhlinin et al. 2009a; Sun 
et al. 2009; Eckmiller et al. 2011) that have again been 
corrected for self-similar evolution (see Figure 8). Khos- 
roshahi et al. (2007) suggest that their sample of FGs and 
"overluminous elliptical galaxies" (OLEGs), also shown 
in Figure 8, are hotter for a given mass compared to nor- 
mal systems, an effect that increases toward lower mass. 
They tentatively attribute this to a lack of cool cores in 
the FGs driving the emission-weighted temperature to 
higher values compared to normal systems. We also find 
a systematically higher Tx for a given mass in our sam- 
ple; however, the factor of ~ 2 systematic error in the 



mass is more than adequate to explain this offset, and as 
shown in the Lx~M relation, it is more likely we have un- 
derestimated the masses than overestimated them. We 
note that Vikhlinin et al. (2009a) and Sun et al. (2009) 
exclude the core of the X-ray emission within 0.15 rsoo in 
their analyses, so cool cores (or lack thereof) cannot di- 
rectly cause the temperature differences seen in Figure 8. 
Our derived M-Tx for FGs is stil consistent with that of 
normal groups within the relatively large statistical and 
systematic errors. 

4.3. Notes on Individual Systems 

JO 133- 1026— The X-ray IGM of this system is the 
faintest in our sample, and we are unable to constrain its 
temperature or mass. In contrast, the BCG contains a 
bright apparent point source with a flux about one-third 
that of the thermal group emission. This source is well- 
fit by a power law with Y = 1.8+ as described above, 
and we conclude that it is an AGN in the BCG. This con- 
clusion is supported by two FIRST radio lobes projected 
8" (17 kpc) and 11" (23 kpc) from the BCG center, with 
a total flux density of 370 mJy. While very bright, these 
sources were missed in our filtering discussed in Section 2, 
which only eliminated FIRST sources within 3". The op- 
tical spectrum lacks emission lines indicative of an AGN 
(Mauch & Sadler 2007). The extended X-ray emission 
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does not peak at the BCG but appears as a clumpy ring 
out to a few hundred kpc. The exclusion of the central 
point source is unlikely to cause this structure, since the 
region we excluded is quite small compared to the extent 
of the diffuse emission. Deeper observations will clarify 
the morphology and temperature of the IGM. 
J 08 15+ 3959, J 1039+ 3947, J141 1+5736— These three 
systems share some morphological traits. The diffuse X- 
ray emission is not centrally concentrated in these shal- 
low observations, and each BCG is separated from the 
peak of the diffuse X-ray emission by more than the un- 
certainty in the center of the f3 model fit (typically 6"). 
J0815+3959 is devoid of an obvious BCG X-ray source, 
while J141 1+5736 has an extended BCG source that pro- 
duces 23% of the flux of the IGM. This latter source 
might simply be the peak of the IGM, and its exclusion 
could explain the unusual morphology seen in Figure 4. 
J1039+3947 has a point-like BCG X-ray source, as well 
as a 17 mJy FIRST source projected 3.2" (6 kpc) away, 
just missing our radio source filtering. J1039+3947 fails 
the Ai > 2 mag criterion using the r2oo value derived 
from the temperature, r2oo,fcT = 916-1273 kpc (see Ta- 
ble 5). It meets this criterion for the smaller value derived 
from the /3-model fit; r2oo,/3 = 809 kpc. 
J0906+0301 — As mentioned in Section 4.1, this target 
was not detected in our 10 ksec Chandra observation, 
yet 25 galaxies around the BCG exhibit a similar red- 



shift with a velocity dispersion of a = 506 ± 72 km s _1 
(Proctor et al. 2011). From the scaling relation presented 
by those authors, we infer Lx ~ 3-30 xlO 42 erg s -1 ; this 
is mildly inconsistent with our 3a upper limit of 4.8 x 10 42 
erg s _1 . Deeper observations of this target will be illu- 
minating, as it could be both optically and X-ray faint 
for its inferred mass. 

J 0856+ 0553, J 1007+ 3800, J14 1 0+4 14$ — Three of the 
four most X-ray luminous FGs, these systems have well- 
constrained temperatures and /3-model surface bright- 
ness profiles. The X-ray isophotes are regular for each 
FG, with the optical BCG projected near the peak of the 
diffuse X-ray emission, although J1007+3800 has some- 
what compressed isophotes on the east side compared to 
the west. J0856+0553 contains no detected BCG X-ray 
source, while the other two have extended BCG X-ray 
sources producing a small fraction of the total X-ray flux. 
J0856+0553 fails the > 2 mag criterion at the upper 
range of allowed r2oo, but is otherwise consistent with 
the FG definition. 

J 10 17+ 01 56, J 1133+ 5920, Jl 153+ 6753— These systems 
appear relaxed and quite compact, with the detectable 
diffuse X-ray emission all within a radius of 190 kpc. 
All three have X-ray centers within 6" (12 kpc) of 
the BCG location, which is within the typical uncer- 
tainty for the center of the f3 model surface brightness 
fit. Jl 153+6753 clearly has compact emission associated 
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with the BCG; we note that this target falls outside of 
the FIRST survey field, but there are no other cataloged 
radio sources consistent with this position. J1017+0156 
has a more extended X-ray source associated with the 
BCG. Jl 133+5920 has no obvious source in the center 
aside from the extended IGM emission, although there 
are FIRST radio lobes projected 3"-4" (7-9 kpc) from 
the BCG, with a total flux density of 11 mJy. With 
the updated r2oo values shown in Table 5, J1017+0156 
strictly fails the > 2 mag criterion, with = 1.8. 
Nevertheless, it is clearly an optically under luminous sys- 
tem (Proctor et al. 2011). 

J 1045+ 0420 — This is the most X-ray luminous FG in 
the sample and also the most distant. The morphology of 
the diffuse emission is strikingly irregular; the peak X-ray 
surface brightness is projected 1.2' (190 kpc) northwest 
of the BCG, and a fairly uniform bright region extends 
from this point through the BCG to nearly rsoo in the 
south. A simple f3 model centered on the BCG or X-ray 
centroid is clearly not a representative surface brightness 
model, therefore the estimated mass and extrapolation 
of the luminosity to r$oo are likely incorrect. There are 
two bright radio lobes projected ~ 15" (40 kpc) east and 
- 25" (67 kpc) west of the BCG, with a total FIRST 
flux density of 160 mJy. The irregular morphology and 
bright radio lobes could be evidence of a recent merger 
or other activity. Deeper X-ray observations will help 
clarify this puzzling system. 

J1136+0713 — This group has a bright red-sequence 
member 0.5 i mag fainter than the BCG and projected 
about 3.8' (440 kpc) to the west. Thus the A* > 2 mag 
criterion is not met at the upper range of allowed r2oo, 
although it is consistent at smaller r2oo derived from the 
/3-model extrapolation. Deeper X-ray observation will 
better constrain r2oo an d determine whether this is a 
real FG. The diffuse X-ray emission is fairly regular and 
peaks at the location of the BCG. No separate BCG X- 
ray source is detected in this target. 

5. SUMMARY 

Fossil groups present a puzzle to current theories of 
structure formation. Despite the low number of bright 
galaxies, their high velocity dispersions (e.g., Proctor 
et al. 2011) and high gas temperatures seem to indi- 
cate cluster-like gravitational potential wells. There have 
been very few FGs with good quality X-ray data observed 
until recently, and their idiosyncratic characteristics may 
contribute to enhance their apparent contradictions. We 
have embarked on a project to assemble a large sample of 
optically identified FGs with a view toward dramatically 
increasing the number of such systems with high quality 
X-ray data. 

The principal observational results in this work are: 

(1) New Chandra X-ray detections were made for 12 new 
FGs, from a sample of 15 optically selected groups from 
the maxBCG cluster catalog with richness ranging from 
9 < iV 2 oo < 25. 

(2) The new X-ray data yielded temperatures for 11 of 
the FGs, ranging from 1.3 to 2.7 keV. From these temper- 
atures and an analysis of the surface brightness profiles, 
we have estimated rsoo ranging from 440 to 710 kpc and 
masses ranging from M500 = 0.3 to l.OxlO 14 M . These 
values of masses and scaled radii are typical for groups 
and clusters in this temperature range. 



(3) The Lx~Tx relation for these new FGs does not devi- 
ate significantly from the expectation for normal systems 
intermediate between clusters and groups, although they 
tend to be more similar to galaxy clusters. 

(4) The Lx~M and M-Tx relations suggest that the 
FGs are on average hotter and more luminous than nor- 
mal systems, similar to the results of Khosroshahi et al. 
(2007). However, the systematic error from luminosity 
correction and mass extrapolation are large enough to 
explain these differences. 

(5) A small number (10-40%) of the detected groups are 
morphologically irregular, possibly due to past mergers, 
interaction of the IGM with a central AGN, or superposi- 
tion of multiple massive halos. Two- thirds of the X-ray- 
detected FGs exhibit X-ray emission associated with the 
central BCG, although we are unable with the current 
data to distinguish between AGN and extended thermal 
galaxy emission (e.g., embedded galactic coronae, Sun 
et al. 2007). 

We conclude from these results that the selection cri- 
teria devised in Section 2 were successful in finding real 
fossil groups, and we have greatly increased the number 
of known fossil groups, a crucial step for further statisti- 
cal analysis. The results obtained from further detailed 
studies of this sample will have implications for current 
and future cluster population studies, and also for cos- 
mology using galaxy clusters. The presence of a popu- 
lation of intermediate mass clusters with small numbers 
of galaxies may bias determinations of the mass func- 
tion which measure richness by galaxy counts. This po- 
tentially biased mass function, when used to set strong 
constraints on power spectrum normalization and Q m 
(e.g., Allen et al. 2011, and references therein), may in 
turn bias these results. Furthermore, such bias would 
also affect the measurement of the mass function red- 
shift evolution, which is used to constrain the equation of 
state of dark energy (Vikhlinin et al. 2009b). In contrast 
to X-ray surveys, where good proxies are being refined 
very quickly (e.g., Kravtsov et al. 2006), the largest cur- 
rent and near- future cosmological surveys (e.g., DES 15 , 
BOSS 16 , J-PAS; Benitez et al. 2009) estimate the mass 
of clusters through optical mass proxies, using some type 
of richness indicator (e.g., Rozo et al. 2009, 2010; Rykoff 
et al. 2011). This dependence on a possibly biased mass 
proxy drives the need to determine the mass range of 
the effect, the relative abundance for this massive but 
optically poor cluster population, and proper statistical 
correction methods for cosmology. 

With these data in hand, we have begun a system- 
atic study of this sample. The addition of optical radial 
velocities obtained for hundreds of galaxies in these new 
fossil groups will allow us to study the scaling relations of 
fossil groups presented in our companion paper (Proctor 
et al. 2011). With our planned deep X-ray follow-up with 
XMM-Newton, we will better constrain the IGM tem- 
perature, luminosity, metal abundance, and halo mass, 
and thoroughly explore the morphology of the hot gas. 
Finally, with additional proposed Chandra snapshot ob- 
servations, we will extend our sample to fainter Lbcg, 
further testing the validity of our selection method and 
pushing the scaling relations to lower masses. 

15 http : //www . darkenergysurvey . org 

16 http://cosmology.lbl.gov/B0SS 
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